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SUMMARY 












A theoretical investigation is made of a ciroular jet issuing 
from an orifice into an air stream flowing parallel to the sot athe 
and in the same direction as the jet. The analysis is divided into two 
portions: 

1. The region immediately downstream of the orifice in 

which exists the potential cone af uniform velocity. 

Ze The region downstream of the end of the potential cone 

in which the jet is fully developed, 

at which the potential cone disappears civing a complete solution for 
the whole jet. 

Numerical mapeenines are made for several ratios of free stream 
velocity to jet velocity to determine the spread of the jet, velocity or 
the axis of the jet, and spread of the surface on which the velociey i 
one half the sum of the jet velocity and free stream velocity. 

Sguire and Trouncer's results are plotted es a comparison with those 


obtained in this analysis, 


Il. INTRKOUUCTION 


the purpose of this analysis is to provide a comperison of the 
characteristics of a round turbulent jet discharging into a moving 
air stream with the known characteristics of a jet discherging into 
air at rest. oth jet and free stream are assumed to have the same 
density and the flow is assumed to be incompressible. 

4% the present time there exists no sstisfactory theories for 
turbulent mixing. However, the interrated forms of the equations of 
motion, along with some dimensional reasoning, can be used to give 
some useful results, 7 

‘The solution of the round jet 4n a general stream has been 
already oarried out by Squire and Trouncer, In their solution the shear 
stress is based on Prandtl's momentum transfer tne oe y amtcl is known 
to be in error, it therefore seems worth while to investigate this 
problem using the integral equations and also making use of experimental 
data for obtaining the shear stress in the mixing rerions,. 

Use of the integral form of the equations was first proposed by 
von ieee; one of the earliest applications was by ee 0 
investigeting laminar boundary layers in a pressure gradient. soreen © 
originally used the intergrated form of the mechanical enerzy equation 
in laminar boundary layer investirations, Liepmann and natitaf Hew 
also applied this method to some simple two dimensional turbulent 
flows, 

In the integrated form of the equations of motion the only 
assumption neccssary is that for the velocity profile across the 


mixing region. This can be made quite accurately anc for e closer 


comparison to the work of Squire and Trouncer the seme velocity prof 


that tney use will be assumed in this solution, 







Ir transformine the momentum equation to the interratea form 
the term conteininr the shear drops out. It is therefore necessary to 
introduce the mechanical enerry equation. ‘This equation, while not 
independent of the momentum equation, is necessary in order to introd 


a term containing the shear stress, 


The following notation will ve used, 


r 
“p 
Mixing Region 
b 
ve  —— Yo 
Jet Axis 
— 
Yb 
Nozzle Exit End of 
Potential Cone 
Fige 1 
x Distance from nozzle exit alonz jet axis 
r Distance from jet axis 
Xo Distance from nozzle exit to end of the potential cone 


a Radius of the potential cone at any point 

b Radius of the outer jet boundary at any point 
b Uiameter of nozzle exit 

u Velocity parallel to jet axis at any point 

a vet velocity at nozzle exit 

u, Free stream velocity 

Ue Velocity on the jet axis at x =x, 

Vp Velocity of the stream normal to the jet axis 


QO Stream censity 


UT Shear stress 


LA 


II, ANALYSIS 
fe Kegion between the exit and the end of the potentiel cone, 


the three fundamental equations are 


Continuity 
roltor' aa SLOP 10) =O 
OX gi 
“omentum 
OU OU 
CM ax TOU SE = par re) 


sechanical energy 


Qu ok +ouv Sb = 4 Berry 


These three equations are intesrated between the 
coundaries of the mixing region and are then 


——— 

| urdr +U,b + Ua Fo —u,b 92 = 0 
QO 

xonentum 


a[u rar +u Wb + Us a¢s - ub 22 =0 


Mechanical enerrcy 


U 
1 nie urdr +4uzu,b +4uda $e ~pudb ge =— | UP Se 
2 aX a 


A funetion is assumed for the velocity profile 


across the jet such that 


u =u, + 4554211 ~ cos. | 


ae a7 Ue 
U=U,t=3 > £(2) 


where 


f-c3 
r = b-(b-a)v 
dr = + -a) dr 


Substituting into the continuity oquation 


$| fuse sg Fyfo-to-abdp |[-6-a)y] 





- +¥,0+ U0 ga 7 u,b a =0 
Expanding | 


S| ulab-b[dy +ua-vi] dre + sab 2) [Feqnay 


+ Ae“ (a- oF flay | nb + usa $y - U.b ae = 
Interrating 
é -ulab-b) - Sla-b)'+ 454° (ab-bYI, 


+ Us (a-b) I] | + Ub +U, a a6 — Usb ge =O 


a = ff (m) ay 
E = { etopdy 





Difrerentiatin. and collecting terms 


[Aa+Bb] 34 a0 +[Ba+Cb] 22 +U,0 =0 (1) 


where 


A= (Ue-UI,+1) 
B= (Ua-Us - I,) 
CG = (Ue- TR lense 


Substitutinz into the momentum equation 


fo [uss Ser Se f)] | @b-b') bee dx 
+U,U,0 +uza oe -uzb 9 
kxpanding and intecratinz 
a [-urdab-by) ~-42(a-0) +u,(U.-U J@b-b) I, 
veya, S520 a, + sgt ol 


+ Up, 0 +UZ0 da. ~ ao =O. 


where : | Fra dy 
I, = [xP cody 


Lifferentiating and collecting terms 


[Da+Eb| 44 *[Ea+Fel oe + ULV, D =O (2) 


dy 


where 

De us =u + £8 (Us-Us)° +21, U,CUe- U4) 

E = (33 - $8) (uu +(1,-21,)U Uae) 
_— (fs 2 As) (Us-Us) +(21,-2T)U, (Ua-Uy) 


The shearing stress now appears in the mechanical 
energy —- and an assumption will be made for the 
shear profile, The shear stress is a variable and could 
be determined from the velocity profile . However, a 
shear profile derived from an already assumed velocity 
profile would heave « greater degree of error due to 
amall errors in the assumed velocity profile. For 
simplicity an average shear profile will be assumed 
_ for the whole region. 

From experiments with a jet discharging into still 
air sheer profiles have beer found as shown in Fig. 2. 

At the jet exit the shear profile is symmetrical. Past 
the end of the potential cone the maxincum shear occurs 

at about % =-067 to .75 and the shear profile has a steep 
Slope on the high velocity side. <An average sheer profile 
for the potential cone region is shown in Fig, 2 and is 


given by 
Y =e UU’ 
which is assumed to be of the form 


Te e (la -Ur) OY) 


wners 


g (7) = 
2 Ts 
In perticular it will ove assumed here that the 


function e(y) can be represented by 






9 (v) =ky?cos* T(R- 
Substitutisg into the mechanical enerzy equation 


id. } [use ust f@){ @b-b') +(a-b)| dy 


a type pee, — yee 


c Ghat 
Le. 
ae 
rc. ae 
ff’ 
Sault. 
\ 

y ia | 
Ty 5 
| Ta 
| , 
| 


ough | (ua-U9@) b- (6-a)@]|-6-a)d4y| 


Expanding and integrating 


ral 5 ur? = 9 Ub Re ath a oa Uy (War U,) I, 


+% Ge u,) I s|[ab- b| 4 +|- bu, + 5 Ug (Us -Ue) Le 
3 2 ! : : 

+ = Ub (Ua- Uy) Le — g (Us-Uy) T,][a-b} 
+uUusyuUu,b + udca go — Usb a. 


= TT (@.-ug Gee ee (Us-U,) 1, Ol 


where 


Je= 4 
I,= vf) 
l= f@g@) 


I,= .7IQI®@ 
Differentiating and collecting terms 
[GarHb|s +[Ha+Jb| ge +uryb=Ka+Lb (3) 


where 


G = Ug -Ug + BUg (Ue-Ub)T, +> U(ue-Up) Ly 
+ +(u.-u,) I, 

H= uf(u.-u.)(E1,-31,) + uy(ua- Us) (a I,-3 1,) 

| + (Us-U) (CS1,- 41) 

T= Ud Ua-Us(S1,-31) + Uy (Ua-us)*(31,-3 1) 
+ (Us-Ur) C41¢ - 415) 

K = 7 (Ua-Us)” [, 

b= er Ca-u,5.(1,-1,) 


Eliminating vy in equation (1) and equation (3) 
[Ga +Hb] $44[Ha + Jb] $b - uz [4a +Bb] $4 
-ui [Ba+Cb| 2 = Ka+Lb 





le-uzA] aa + JH-usB]/o Sz + oG8-| 
J-urC] 42 =Katlo 


J, [cat +2Gab + ob] = Ka+Lb 


where 


of = SUA 


_ H-uzB iene 
2 j G = il ZZ 


Ze 
Using the chain rule of differentiating, a0 = ah gb 


the above equation becomes 
l2xa da +2Ga +2Gb 44 +2rb] ao = Ka+ Lb 


(Ka +Lb) $¥ = (2«a +2@b) 42 + 2 Raggee cae (4) 


Eliminating v,, in equations (1) and (2) 


Da+Eb] 3s Ea +Fb] die _ 
fe) 


b 
-wilads + [e-ualeage +e da] deus «c 
Inte rrating 


507+ 2e€ab+¢b'+Z =0 (5) 
where 
= Sauna = Sie Oss _ F-U C 
: oS) 0S - —— 


Z = constant of intecration 





Differentiating with respect to b 


26a 54 +2E€Q+2Eb 2% nC =O 


aa _  €42+Ob6 (6) 
db Eb + 5a 
Substituting equation (6) into equation (4) 


Ka+Lb| 2% aX = aa +23 QC EA+ > garb). 2 (30 +270 


Caer a 


_ (8B -a€)a+ (7S - «dab +(wE-oA) b- 








sei _ 

2db Eb + 6a 
1 ax 7 A2b Ask 

ome > A+ acer + ast t 


= 
= 
@ 
4 
o 


i §G- «Ee 
| 


&K 


n _ (&8-<*)(6B- x€) 
2 §Cek — SL) 


X= KE(LE-Ke)+Ko(KG-Led) + (ae - $8) 
: K(L& -Ké) 


From equation (5) 


Sat+€b = -Z48 + bi(e*-o 5) 


td dy _ —— ALb Aa 

fab ea 

Integrating 

x a AV ~-28 + CE*-45) —_bdb- 
=)\ b- eee As | bab 


This expression can be integrated and Zzives 


oi hee a V-Z8 V~28 +b (E*-$6)_ 


— Oe- eS 
: —s bYe,+! -¢ - CAG 77) 
+ Azle, a = Ge —C 2) == lo Og bVc.et + C on +10(C, Sara 


"3(C+l) 





where 
_ § 
Co= TE-KE 
_ _k*z68 
@s wey" 
@ RG Se 
2 (L&E-Ke)y* 


Y2constant of intercration 


Ga 


ty 


Y can be found by the boundary condition that x =0 
where b= d/>. Obtaining the parameters x and »b in 


dimensionless form the result becomes 


- 4h - &) ($8 - €*) 


x 

a= 2A, (e- ) - eg 

+ ay ~ 432— (08 -@ 
ps -eE* 


+Aateb fa Gey goer ng GO NE MT ca] 
Cy+ | B Be Se a) afc, +1 09 acne * ra aan Ze\ cI} oe 


A3Co | - ia, og C.+! rhe, +D1CE)~c ] +2E4 
az ‘ae A Ne. a; + Morar aT re 


The constant of integration 2 is obtained from the 








boundary condition that a=b=Y at x=0 so 


Z = -£u,(U-Us) 


‘The only remaining constant to be determined is Ae 
which appears in the shear term, From experimental data 


©) ss 4s found that near the 


obtained by Liepmann and Laufer 
jet exit GO ee = 0.0072, From experimental data 

obtained by odbesin’ a€ is Péurid that in the fully developed 
region GM) eax. = 0,022, Both of these results are for 


the case of a jet discharging into still air or as =O 
Qa 


Differentiating our assumed function it is found that 


_ 
yl 
d 7 . pes > = = a=, 
ae a8 FLY PPA HO 
iy 3 -— Ce oo ol a2 7 
sf dip i, , FF > — ; 
e .& é ° ; 


The constant must be given some average value such that 


G(7) max occurs at the 





point where 7=0.692 


4 nan f oF the average sheer profile will be between 0.0072 
and 0.224, 

From experimental data it is known that the end of the 
potential cone is at a point about 9 radii downstream of the 
exit of the jet for the limiting case where Gt =0. in order 


ee 
for Y, = 9 when Bone and a=0 
Fh (Pgy = 0.01787 


and fe = 0,078&8 





8OUG 


: 
Lo 


B. Region downstreem of the end of the potential cone. 


At the point downstream of the jet exit where x =x. 
the potential cone disappears and the variable, a, no | 
longer appears. Downstream of this point the velooity on 
the jet axis, u, , will become a variable decreasing with x, 
The assumption for the velocity profile in the revion of 


the potential oone was 
Ua-~U b-r 
Cee iA + Barve || - COS 7 (8-5 )] 


a-=0O and UL= Ua 





At X= Xo 


ree Cle + Ke=Ue| | Cg (PS) 


This reduces to 


u = Uy+ Page |i + cos T | 


or 
(A. -U 
Ue Uy + St $y) 
where 








The continuity equation in this rerpion then becomes 
d oan UWe-~Up > ) f d | b b db ee 
d YX Uy D |p dpa i lo b {. (mM) AM. a3 Uy =o 1b EB ay = U 


This can be reduced to 


1,V $+ 21,6 gt 45. ae 


where 


= [toogp 


I, _ 
V = U.-U, 


Substituting in the momentum equation 
flush pap * Us aeuyb’] tow dun 
+a (Uc-Ub) re pf “w)du.| +uyy,b- ~U ob 38 =O 
_ Which reduces to 


41ju,V+1V4 $e P *[2Tou + V]JbS+ 2uy vn =0 


where 


I =| uf qu) dw 


_ 





In the mechanical energy equation the sheering stress is 
given by 

v= O(Ue- Uy) GG) 
The funetion 9 (pr) must be such that the sheer profile will 


resemble that illustrated in Fig. 2. This will be 
L je 
gy = 2 P| a fwdy 
? = constant 


Substituting into the mechanical energy equation 


£ luz By dys +S Up (Ue- Ae Panty Ap 





+2 Ue (UU bb, pf Gu)dye + Uc SOT uf aod] 
+ UU. b-uzb go =—b Ue-Uey | q(u) + QW) A 


Which reduces to 


| | 31 ,ueV a5 $1,U.V + = I, v> | gb 
du. 
+[21.utb +31 usbV +21, bV7] SS 


+ Cran ore = ~Il,V° (10) 
where 


le =f pf wd 
Le =[ maggot qu dp 









Eliminating v. 


» in equation (8) and equation (9) 


4.1,uV +Lv 1S % l2 Tou, a IVb wv 
- 2u, LIV 92-4 15b 3 = 0 


db ae Tou, +1..V dV 


es =a ee 


bh 2Touy + TeV V 


Inte grating 


log b =— + log|2T,u,+1,,V{- . log + log G 





— | } 


Aly (2 TU. +1,V) 


db dew ae ToUetlioV : av. 
dy , V Clout Te. V 
db a C(L5u.+1,V) AV 





2 


dy. V (2 guy +1,,V)* d¥ 


== 
awe 


Elimineting v, in equations (8) and (10) and substituting 
for ae above we -et a single equation which is 


Toluct+loluuVralvo dv _ 


Ly 
fr —_— zon V ° 
(21 uy + TeV) JX 


$ 
This may be integrated to sive 


1 oe 


eee aap) 
SC, V 





-2 (2) ¥ : -@) 


The oonstant of integration C, may be found by the 


boundery condition that at x=x_, V = UU, 


eo 
oF 


¥, \) 

s = q, +2((3)) - me gi 
t -(% | “~%  g Inte 

“LI, r) ~~» © 


LA 
= cr, pe *ee 
where 


“u ~§-2 8.) u, 
. Ua—U, 


The constant of interration C, may be evaluated from 





these same boundary conditions and from equation (5) when 
a=0 


Substituting these values of C, and C, the finel 


solution is obtained. 


Kate (ey- as || 2 é1, A + T(1- 4e)| 
He Ver ree bt 2) 


(a 


3/2 o/s. 7. Ya _ S/2 
ae RE: | Q~ | + |Saks G La 
. & Lee $b Gls1,, 


x] 2 (gage) + Ge" ge 7) ‘q 


where 








= VJ 
$0 
Equation (11) becomes 
bY ee Sr) ed ote 
(ay Ea -\\f2i, + T,, (‘4s Me -) E+ ») Ae $+1-I,~ a 


x, may be found by solving for b when @=Q in 
equation (5). Using this value of b find the corresponding 
value of x in equation (7) which will be x,. 

Equation (13) cives x as a function of (u.-u,). 


Equation (14) gives »b as a function of (u-u,). It does 


not seem necessary to cet an explicit expression for x 
and »b. By assuming values of (u,.-u,) correspondinr values 
of x and »b are found and b is plottec against x in 
Figs Se 

Now the oonstant l appearing in the shear stress 
expression must be evaluated. In the fully developed region 
downstream of the end of the potential cone J (pUma,= 0.022, 
From the assumed function of 4()it is found thet 


BWWrax = 0.4382 or ¢ = 0,05116, 


III, DISCUSSION 


The results plotted for the region containing the potential 
sone ere the results obtained by Squire and Trouncer, In using the 
derived equation for x and b (equation 7) in the potential cone 
recion the final results obteined did not appesr physically reason= 


able, particularly at higher velocity ratios of free stream to jet. 







Due to the enormous amount of computation necessary to recheck these 
results Squire and Trouncer's results were used in this recion, 

The very close comparison of the two analyses in the fully 
deve loped veahen vuaill seem to indicate that the integrated equation 
solution that has been used here should give feir erreement in the 
potential cone region. 

The basic anelysis of the potential cone refion is known to be 
oorrect through equation (5). This equetion zives ea relation between 
the spread of the jet and the width of the potential cone, fhe 
remainder of the analysis within this same region or the computatélien 


may be in error, However, the method of solution igs outlined for any 


future investicetion of this problem. 
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g2reeeiy I 


Evaluation of dimensionless integrals appearing in the 


equations of motion, 


The dimensionless interrals appearing in the equetions 
of motion all may be expanded into integrals of forrs found 
in standard intogral tables, Following are the results of 


ie J fopdy : \fi- cas Tr | diy 
-- J rekeady : [ [i-cos tl pdy 
a ) F@dy =| [I-cos ty] dy = — 1.500 

Ly J gf (pay -| [i-cos ty] pdy= - 15s 

LJ Po) =| [i-cos ty] dy = -2.500 

I. aia vody={ fI-eos ty] dy = - 2.014 

i | geo! mah= | feyy’cos*rr(y-4) sm ny dy =7.068TSxi0 Ye 


-|.00 


I 


~ OFFS 


a =f qootaddy-[Rytcaste(y-4)si tydy =- .04305xI0 7 


bee Fd. -[(i+costuludu = 0.298 

L.-['Btoadu ={ tu cosmu] wd = 0.346 

i | uf Gdn ant +COS reps] an Gir = 0.485 
Le-{ g (AL) Mu oyu 


=|nl sin mal + cos trya|| 4 as de (Cos TA 1A Stn Thy = 


=O) {O) Sy Sisa] 


APPENDIX II 
Integrations necessary for solution of the equations, 


A. Integral appearing in the equatior for the region 
containing the potential cone, 


[=| b db 


Katbb 
Substitutins for a from e.uation (5) 


1-{ b db 
ae? @ eG a Se ttb 


which reduces to 


= 


F C0) b db 
“Sy 4+ 


where 
(., <a 
° ~6hL& -KE 
CC . —— 
7 (LS -Ke)* 
C = K "(> 8 =e) 
2 (Ea— Kee 


Rationalizing the denominator 


(b* bvfo—c.b* ) db 


= © 
I : or (Cee! arene 


L= co] oi a eS a Ed 


b (@.+1)-¢2 lb” (C +1) — E> 
2 b db 
> bt (Cc, + -C.t 
lb db 





C, 
p Pg 2D | I 
am C4 ovVC.4 ~C, bVc,+! aaa 





Tecra ceil 
JES" bY GH 6, VGH +6, 2G.) 





= xt | 20 Tea pose lag (oes? ~c,) - one 





7 C, eames | —C 
~ 2@+ [26 * Year !09 yvecr —_ 


Afci-c.bt » ab 
b(c,+1)-¢? 
{ et c*-¢,b = c.R: 
-2¢, bdb = 2C,RodRo 


bdb = me, & Re 


7 | VC. Ro (- RedRa) 
(C.+0(S 9), eee 


Vcx Re dRo 








(c, «(ieee 
C 
= -7e 2 ee dR, 
RAC + = Jak Rov Cat! + 


= ler. es hog (ReVGeT = Se) _ a (RNG Hye 





log (c,ti)(c,*-c,b) = : 


— 
Se Nl a 
2(Cy +1 (2 N ; 


B, Integral appearing in the equation for the region 


downstream of the end of the potential cone. 
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An investigation of 
the round Jet in a moy- 
ing air Stream. 








